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Abstract 

We consider the quantum spin hall state in a three orbital model due to certain loop current 
order induced by spin-dependent interactions. This type of order is motivated by the loop current 
model which is proposed to describe the pseudogap phase of cuprates. It is shown that this model 
has nontrivial Chern parity by directly counting the zeros of the Pfaffian of time reversal operator. 
By connecting to the second order Chern number, we explicitly show the singularities of wave- 
functions and how they depend on gauge choices. In this case, it is shown that the Berry phase 
can be mapped to Nonabelian instanton. 

PACS numbers: 
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I. INTRODUCTION 



In recent years, there are a lot of interests arise in the topological property of band struc- 
tures, such as quantum anomalous Hall effects (QAH) and quantum spin Hall (QSH). All 
these topological nontrivial states are inspired by the integer quantum hall effects, in which 
it is shown that the quantized hall conductance is also the Chern number of underlying band 
structure 111. Haldane introduced a complex next nearest neighbor hopping in a 2D lattice 

This leads to an extensive study on QAH [3|, It was also shown that Chern number is 
the only topological invariant which characterize time reversal breaking Hamiltonian 



In our previous paper 



5^. 



6|, it is noted that Haldane's model can be realized in the lattice 
model with some special time reversal breaking loop current ordered states due to the 
nearest neighbor interaction. The loop current ordered states without changing translational 
symmetry were originally proposed to arise as broken symmetry states due to interactions 



in a three-orbital mode 
famihes of cuprates 
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for underdoped cuprates [7H9| and have been discovered in several 



-Il2|. These observed loop-current states, however, do not lead to the 



quantized anomalous Hall effect (QAHE) or nonzero Chern number because they preserve 
the product of time reversal and spatial inversion symmetry. But other possible type of 
loop current order will lead to QAH. Also in [S], the discussion of the singularities of wave 
function helps to reveal the topological obstruction in QAH states. 

In this paper, we generalize these results to include spin degree of freedom and consider the 
ordered spin loop current states due to the nearest neighbor interaction. This will restore the 
time reversal symmetry and also make the Chern number trivially zero. Although the first 
order Chern number is zero, it is pointed out that the time reversal invariant Hamiltonian 
is characterized by the second order Chern number [3]. Many yeas later, in the seminal 
paper [l^, Kane and Mele introduced the Z2 invariant to classify the two dimensional time 
reversal invariant model, which also lead to the suggestion of 2D and 3D topological insulator 
[15H17[. These two classification of time reversal invariant model seem to be contradicted 
to each other. This contradiction can be resolved by the dimension reduction method as 
in 18|. Another more geometric method to connect the Chern parity to the second Chern 
number can be found i n | l9l| . This argument is very similar to the one used in the Witten's 
SU{2) global anomaly |20|. We will first follow Kane and Mele to construct the Z2 invariant 
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directly. After connecting the Z2 invariant to the second Chern number, we will analyze the 
singular points of the wave function of a two band model example in detail just as in the 
QAH case 6|]. For the two band time reversal breaking Hamiltonian, it is well know that the 
f/(l) Berry phase configuration looks like a monopole solution {21]. Similarly, for the two 
band time reversal invariant Hamiltonian, we found the SU(2) Berry phase configuration 
can be mapped to the self-adjoint instanton solution of Non-abelian gauge theory after some 
suitable coordinate transformations. At last, the three orbital Cu-0 lattice model will be 
discussed. Again we will directly construct the Z2 invariant to verify it is a QSH state. 



II. SPIN LOOP CURRENT STATES IN THE COPPER-OXYGEN LATTICE 
MODEL 

We consider the two-dimensional lattice with the structure of the copper oxides lattice 
as shown in Fig. ([T]). There are three orbitals per unit-cell, which are the (i-orbital on the 
copper atom and the px and Py orbitals on the oxygens. The minimal kinetic energy operator 
with a choice of gauge such that the d orbital is purely real and the p^ and the Py orbitals 
purely imaginary is 

Hke = it dl{sxPx,]<i + SyPy^k) - t' SxSyP^^^]j)y,\, + h.c. (1) 

with Sx = sm{kx/2) and Sy = sm{ky/2) for a lattice constant taken to be 1. 

In the loop current model of pseudogap phase, one usually only consider the charge 
channel, thus is equivalent to consider spinless fermions. To consider the spin Hall effects, 
we also need to decompose the interaction term in the spin channel. As we know from 
the previous paper, the only topological nontrivial state is obtained by considering the 
interaction between the p orbitals as 

Hint = ^ ^ VUp^i^uUpj^a'- (2) 

Here i,j labels the lattice sites and a, a' labels the spin. 

In the spinless case, the the above interaction term can be decomposed by using the 
operator identity niUj = — 1(| Jjjp — — nj) with charge current operator Jij = i{clcj — c^jCi). 
In the mean field theory, one of current in the current interaction term can be replace by its 
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Figure 1: Schematic figure of loop current in Cu-0 lattice 
expectation as {V/2){Jij) = ir, thus one finds an interaction induced kinetic energy term 

H-nt = ^rc^Cypl^Py^k + h.c. (3) 

If r 7^ is a stable state, it describes loop currents flowing clockwise (or anti-clockwise) 
around the oxygen's in each unit-cell as shown in Figure ([1]). This is one of the five possible 
loop-curent states with non-overlapping loops in the Cu-0 lattice all of which preserve 
translational symmetry. In Figure ([1]), the flux has one sign in the square formed by the 
nearest neighbor oxygens which surround a Cu and the opposite sign in the square formed 
by the nearest neighbor oxygens which do not surround a Cu. Therefore the total flux is 
zero. 

For the spinful fermions, the above operator identity can be trivially generalized to the 
interaction term which are diagonal in spin indices. Following similar method, the den- 
sity coupling can be rewritten as spin currents coupling by operator identity rii^rjnj^a = 
- n^,a - nj^a) whcrc the spin current is Jij> = «(4<^Cj> - cj^^Q,,^). 

For the off-diagonal terms in spin indices such as ni^^rij^i + rij^^rii^i, the decomposition is 
more complicated. We can begin with the product of the spin current of two different spins. 

= StSj + SrSi - T^T. - Trr; (4) 

Here = 4^0^, S~ = cl^Q^ and S- = ^{rii^^ — rii^i) form the spin SU{2) algebra, and 
= cl^cl^, T[~ = Ci^Ci^ and T/ = \{ni^^ + rii^^ — 1) form the charge SU{2) algebra. Making 
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use of the following identities 

s,-s, = ^{s^sr + s.si) + s^s; (5) 

T. • T, = ^{T^^Tr + Trr;) + T,^T^^ (6) 

we find the following decomposition. 

n,,n„, + n,,n,, = J.,,J,,, - 2(S. ■ S, - T. ■ T,) + + n,) (7) 

Other than the spin current couplings, there are also more complicated terms like spin 
and charge operator interaction between neighboring sites. Since we only want to consider 
current order here, the spin and charge interaction term will be ignored. As before, in the 
mean field level, the current can be replaced by its expectations, then we find a new kinetic 
term 

Hint = iraC^Cypl^^^^Py,^^^ + h.c. (8) 

If the expectation of spin up and spin down current are opposite, this loop current ordered 
system is time reversal invariant again and is possible to give rise a quantum spin hall states. 

We will consider the QSH state of the 3 orbital Cu-0 model and therefore the singularities 
of wave functions of the model with the Hamiltonian H = Hke + H'^^^. Before we do that, 
let us first consider the simpler case of a two band model. 



III. TWO BAND MODEL WITH SPIN 

In [gI, we have discussed a simple two-band Hamiltonian on a square lattice from a wave 
function point of view. 

H = H ■ a, R = (m + cos + cos ky, sin k^, sin ky) (9) 

This Hamiltonian breaks time reversal. For < m < 2 or —2 < m < 0, the Chern number is 
+1 or —1 respectively for the lower band. In the two band model, the Chern number simply 
equals to the winding number of the mapping R/|R| from a 2D torus to a 2D sphere S"^. 
In this case, geometric meaning of Chern number is quite clear. If we change one of the Ri 
to be —Ri, then the orientation of the above mapping is reversed, then Chern number will 
change sign. 
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Now we can introduce spin degree of freedom to enlarge the above Hamiltonian to become 
a 4-band model and restore time reversal symmetry. Then we can show this model has non 
trivial Chern parity or Quantum Spin Hall (QSH) effects. The only time reversal breaking 
term in Eq is -R3. The time reversal is restored by multiply a spin operator, thus we find 
the following model 

H = H ■ a, R = (m + cos + cos ky, sin kx, a + Sz sin ky) (10) 

Here is the z component of spin. For = 1 or spin up, we have Chern number equals to 
1 and for = — 1 or spin down, we have Chern number equals to -1. Here we also insert 
an arbitrary constant a there in order to lift the degeneracy of the spin up and down bands. 
In matrix form, the model is 

H = Ri{ai + R2{a2 /) + 0(^3 ® /) + sin A:^. (as ® s^) (11) 

Here Ci^2,3 and s^^y^z are Pauli matrices and / is 2 by 2 identity matrix, Ri = m+cos fc^+cos ky 
and i?2 = sin ky. 

The spin up and spin down two bands are still degenerate at ky = 0, ivr. So one cannot 
distinguish the Chern number between the spin up band and spin down band. The total 
Chern number is zero. But for certain choice of parameters, there is also Z2 topological 
invariant. There are several different way to get this Z2 invariant. Here we use the approach 



proposed in Kane and Mele's original proposal 1^. In this method, we define the pfaffian 
P(k) = Pf[(Mj(k)|0|Mj(k))]. The zeros of -P(k) always come in pairs, if k is a zero point, so 
do — k. We can count the number of the pairs. Since two pairs can annihilate each other, 
even number of pairs are topological trivial, odd number of pairs is QSH state. 

In our model, we define the pfaffian P(k) = ('U4.(k)|0|-U|(k)). Here u^^i are the eigenstates 
of spin up and spin down for the lower two bands. G = I®{iSy)K is the time reversal operator 
and K is the operator to take complex conjugate. For concreteness, we take m = — 1 and 
\a\ < 1 and define Rs^^ = a + sinky and R^^i = a — sinky and R^^i = {R\ + -R2 + -^3,t,4,)^'^^- 
Since will always change spin up to spin down, we can omit the spin part of wave function 
and only write out the orbital part of wave function explicitly. 
We already know the wave functions, one way to write them is 

I . 1 (r^^-R^\ , . 1 ( Ri-^R2 , 

'""^^ V2i?t(i?t - ^3,t) \ Ri + iR2 } ,/2R^{Ri + Rs,i) [ -R,,^ - i?/^ 



The above wave functions are well defined if Ri and R2 are not equal to zero at the same 
time. 

First we consider a special case Ri = R2 = 0. We will that show that if Ri = R2 = 0, 
then P(k) = 0. Ri = R2 = only happens for = and ky = At these two points, 
Eq.( IT2l) may be ill defined. In this case, it is easy to find the eigenstates if we go back to 
the Hamiltonian. The Hamiltonian for spin up and spin down at these points are actually 
diagonal matrices 

((a + sin ky) \ I (a — sin ky) \ 

= r (13) 
—{a + sinky) I \ —{a — sin ky) I 

At these two point sinky = ±1. For example, for k^ = 0, ky = 7r/2, then we find = (°) 
and 1^4.) = (q), from this clearly P(k) = 0. Similar thing also happens at kx = 0, ky = — 7r/2. 

So far we have found one pair of zeros. Now we want to show there is no other zeros. 
Thus this is the only pair of zeros and the number is odd, so we have QSH state. Since 
we are only looking for zeros, we can drop the normalization factors. For cases other than 
Ri = R2 = 0, Eq. ([T2D is well defined, thus we find the pfaffian as 



{ui\Q\u^) oc {Ri + iR2, -R3,i - Ri) 



Ri - iR2 

R, (i?3t -R3i-Rt-Ri) + tR2 {Rst + Rsi -Rt + Ri) (14) 



Now we show if we assume (m4,|0|mi-) = 0, then it will lead to some contradictions. We 
can distinguish 3 different cases. Note that we always have the following identity 

R]- Rl = Rl,i - Rl,i (15) 

First, if i?i = and R2 7^ 0, then {ui\Q\u^^) = implies R^ — Ri = R^^^ + R^^i- Combining 
with Eq. ( fTSl) . we find R^ = R^^^ and R^ = —Rs,^- But since R2 ,we should have 
> |i?3t| and |i?4.| > |^34.| which is contradict with above equations. 

Similarly, If i?i 7^ and R2 = 0, then (m4,|9|u-|-) = implies R^ + R^ = R^^^^ — R^^i- 
Combining with Eq. (ITSl) . we find R^ = R^^^ and Ri = —^3,4, Again since Ri 7^ 0, we should 
have |i?4-| > I-Rsti and |i?4,| > I-R34.I which is contradict with above equations. 

At last. If both i?i 7^ and i?2 7^ 0, then {ui\Q\u^) = imphes R^ + Ri = R^,^^ — R^^i 
and R^ — Ri = R^^^ + Rs^i, which means that R^ = R^^^ and R^ = —Rs,^- Again since 
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i?i 7^ and -R2 7^ 0, we should have |-Rt| > I-Rsti ^^nd l-Rj,] > I-R34.I which is contradict with 
above equations. Therefore, in all the above 3 cases other than R\ = R2 = 0, we always find 
P(k) ^ 0. Thus ky = 0, = ±7r/2 is the only pair of zeros. 

The Z2 invariant of this two band toy model can be determined by a much easier formula 
proposed by Fu and Kane, provided the constant a is zero. In this case, the model is also 
invariant under the spatial inversion, then the Z2 invariant is the product of the eigenvalues 
of the spatial inversion operator at the 4 time reversal invariant momentum (TRIM) points. 
We assume that the pauli matrices in our model describe pseudo-spins, and the spatial 
inversion will exchange the two sublattices. Under these assumption, the inversion is given 
by the operator P = ^ I- It is easy to verify that H = PHP. The eigenvalues of P is 
sgn(i?i) at the 4 TRIM points, which are one +1 and three —1 or vice versa. Therefore, the 
Z2 invariant is —1 which implies a topological nontrivial QSH states. 



IV. RELATION TO SECOND CHERN NUMBER AND SINGULAR POINTS OF 
WAVE FUNCTION 

The above Z2 invariant is directly constructed from the matrix element of time reversal 
operator 0. But its geometric meaning is not very clear comparing with the Chern number 
of the time reversal breaking Hamiltonian. As pointed out by B. Simon many years ago, the 
topological invariant of time reversal invariant Hamiltonian system is the second order Chern 
number 



13| . The nonzero second order Chern number corresponds to the nontrivial elements 
of 7r3(sp(2)) = 7r3(S'^) = Z. This result seems contradict with the Z2 classification. Actually, 
the second order Chern number only works for a Hamiltonian in 4D. If applying it to 
lower dimensions, one can treat extra spatial coordinates as adiabatic changing parameters. 
Following the similar arguments as Witten's SU{2) global anomaly [l9|, one can show that 
for lower dimensional models the 2nd order Chern number reduces to a Z2 invariant and 
the topological nontrivial one corresponds to odd 2nd order Chern number, which also 
corresponds to the nontrivial elements of 714(5^(2)) = n^^S'^) = Z2. 

In the more general 4D model, the geometric meaning of the topological invariant can 
be understood much easier. In this case, the nonzero Chern number is also the topological 
obstruction which prevent us to define the wave function and Berry phase globally on the 
base manifold. We can understand this in a exact parallel way as in the first order Chern 
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number of QAH state. To this ends, we consider the following model. 
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H = J2R^^^ (16) 

i=l 

with Fj = 0"! (g> Si for i = 1, 2, 3 and r4 = cr2 (g> /, Fs = as ® /. The connection with the 
model of Eq (jH]) discussed above can be obtained by taking 

R = ^ sin kx, sin ki, sin ^2, sin ky, m + cos kx + cos ki + cos k2 + cos ky^ (17) 

and then treat the the momentum fci 2 as adiabatic parameters. 

The energy bands are E = ±R with double degeneracy. Here R = ^JeLi^I To be 
explicit, we only consider the lower two degenerate bands. The wave function are given by 

\^p,) = I (-R3 + tR4, -Ri - tR2, i?5 + R,oY (18) 

2R{R + R5) ^ ' 
|V2) = ^=i^=(-i?i+«i?2, i?3 + ^i?4, 0, R^^R) (19) 

For convenience, we perform the calculations in terms of differential forms. Then the non- 
abelian Berry phase is Aij = {ilji\d\ipj) , which is a su{2) Lie algebra valued one form. Aij is 
traceless, thus can be expanded by Pauli matrices as A = ^21=1 Then we find 

1 .R^dRi — RidRii + R2dR^ — R^dR2 

= ' WrTW) 

,0 R'idRi — RidR':! + RAdRo — RodRi 

r(rJ,) pi) 

,0 R^dRo — RodRi + RA.dR'i — R'^dRd , , 

The Berry curvature is defined as F = dA ^ A t\ A. Here F is also traceless, thus TrF = 
and the first order Chern number J TrF = as required by the time reversal invariance. 
Expanding F by Pauli matrices as F = F°'aa/'2,, we find the following 

F' = -?±3lA' ^A' + ^JdR, A dR, + dR2 A dR,) - . i^ + ^^M'^dR, ^^3) 
R R'^ R 

F' = -?^A' ^A' + ^JdRs A dR, + dR, A dR^) - + ^^f ^ (24) 
R R R 

^ + ^'-A' AA' + ^AdR, A dR2 + dR, A dR,) - + ^^f ^ (25) 



R . X . , . ^2 

The the second order Chern number is related to the winding number of mapping from 
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to S\ 

1 



C2 = 



87r2 
3 



a 

J J2^-'^y^'RidRi ■ ■ ■ dRi-idRi+i ■■■dR^ = -l (26) 



877^ „ 

1=1 

To ease the notations, we omit the A from now on. 

One can see that |'(/'i,2) is ill defined at the points satisfying i?j = for i = 1, ■ ■ ■ , 4 and 
i?5 < 0. To cover the whole manifold we can introduce another wave function with different 
gauge choice. 

I^^i)^ = {R - ^5' 0' -^3 - ^i?4, -Ri - iR2) (27) 

\^P2Y = „ ^ (O, R - i?5, -Ri + ii?2, i?3 - iR^ (28) 

\J2R[R — its) ^ ^ 

1-01,2) s-re ill defined at the points satisfying i?j = for i = 1, • • • , 4 and R^ > 0. Other 

than these points, both of them are well defined and are related to each other by a gauge 

transformation. 

-i I R^ + iRi Ri + iR2 \ 
^0 yR^- iR2 -Rs + iRi J 

Here Rq — ^y Rf + R^ + R^ + Rl- It is easy to see that g is unitary and detg' = 1, thus 
belongs to SU{2). Now the Berry phase and Berry curvature are transformed as — 
g~^Ag + g~^dg and — g~^Fg. The second Chern number can also be understood from g. 
Let IV') and {^p^) be defined on two 4D discs Di and D2. The boundary between Di and D2 

are 3D sphere S^. We know TrF^ are closed, thus can be locally written as TrF^ = du^. Here 
cua = Tr{AF — A'^/3) is Chcrn-Simons 3-form. is not gauge invariant and is transformed 
as a;f = — Ti {g~^dg)^ + da2- Then the second Chcrn number is also given by 

- I lAig'^dgf] = — ^ [ -^y"(-iy-^RidRi---dRi_idRi.i---dRi = -1 

This exphcitly shows that the Chern number is related to 713(5'^). 

To show the relation with the instanton solution of Non-abelian gauge theory, we first 
map the northern half of S'*' to by a stereographic projection as follows 

n = —Ri, Inverse: Ri — , -. R5 — , ^ = (30) 

i?5 ^/WT^ VR^ + 
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for i = 1, ■ ■ ■ , 4. It maps the north pole of to the origin or M'* and maps the equator to 
the infinite boundary of R^. After this mapping, the Berry phase can be written as 

. r^9'^df,g 1 
An = —t , q = -{ri + triCTi) (61) 

which can be verified that do not satisfies the self-adjoint condition F^j = ±^eijkiFl^i. Since 
it has the same Chern number as an instanton, one can always find a continuous coordinate 
transformation to connect these two solutions. Here we can slightly generalize the above 
stereographic projection as 



R, = /(r)r,, R, = ^ R^ - f\r)r\ 2 = 1, ■ ■ ■ , 4 (32) 

Then the resulting Berry phase is 

P(r)g^^dag I, 
■""^"' g + fiVg-f^Mr^ - S = -ir,^^r..i, (33) 

Then it is easy to verify that by taking /(r) = Eq. (ES]) will become the standard 

instanton solution = —i-j^^g^^d^g. This is closely parallel with the fact that the U{1) 
Berry phase of a time reversal breaking two band model is the same as a magnetic monopole 
solution. 



V. THREE ORBITAL COPPER-OXYGEN MODEL WITH SPIN CURRENT OR- 
DER 



We can apply the same method to understand the three orbital copper-oxygen model 
given by Hamiltonian H = Hke + H' as discussed in section [Tll Then the Hamiltonian can 
be written in the matrix form as 



H 



I 



\ 



iR\ —iRi + a 
—iR\ iSzR'i 
iR2 + a —iSzRs 



\ 



(34) 



with Ri = sin-^, R2 = — sin ^ and -R3 = rcos^cos-^-, note that R3 is always positive 
for — TT < kx, ky < IT. Sz is the z component spin. Here we also introduce a small positive 
constant a in order to lift the degeneracy of spin up and spin down bands. With the spin 
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degree of freedom, the Hamiltonian is actually a 6 by 6 matrix as 
/ n „■ D „ I? I ^ \ 

ii?3 
-iRi 

The energy band can be solved from the following cubic equation 



H 



iRi —iR2 + a 
-iRi 
zi?2 + a 



^ 



(35) 



- (a^ + R^)E + 2aRi(s,R3 







(36) 



Here R = ^y Rf + R^ + R^- At general momentum, the last term of above equation is 
different for spin up and spin down, so the degeneracy is lifted at these momentum. Now 
consider the 4 TRIM points (0,0), (0,7r), (vr,0) and (7r,7r). At (0,0), Ri = 0, and at other 
3 points, i?3 = 0. The last term of Eq. (I36l) is always zero at TRIM points, thus the spin 
up and spin down bands degenerate at these points as required. 

As we already know from jol, the top band of the Hamiltonian Eq.f l34p with a = and 
Sz = ±1 has Chern number ±1 respectively. Since adding a constant will not change the 
Chern number, for nonzero a, the spin up and down of the top bands of Eq. fl34l) still have 
the Chern number ±1 respectively. As in the two band case, we want to show that the top 
band has nontrivial Chern parity by counting the number of zero pairs of the Pfaffian -P(k). 
Again, in this case, the Pfaffian is just a number P(k) = (M4^(k)|0|u-f-(k)). 

Suppose the 3 bands are £'1,2,3- It is easy to verify for small a the 3 bands are always 
separated by gaps. We will only consider the top band E^. Since the Hamiltonian is traceless. 
El + E2 + E3 = 0, thus we must have E^ > 0. 

As before, there are many ways to write the eigenstate, we can write it in a particular 
gauge as 



\us,{k)) oc 



aE^-Ri{szR^)-iE^R2 



\ 



\ 



El 



(37) 



i{E^{szR^)-aRi)-RiR2 
Rl 

The wavefunction \us^) is well defined if its 3 components of Ri are not equal to zero at 
the same time. As in the two band example, the zeros usually occur at singular points of the 
wave function. Therefore, we first consider the special momentum point k* where \us^=+i) 
is not well defined or its 3 components are all zero. From the real part of 2nd component, 
we have R1R2 = 0. If we have Ri = 0, then the 3rd component requires E3 = which is 
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impossible, thus we must have R2 = or ky = 0. Suppose we consider the case that -Ri > 
or kx > 0, then we have E3 = Ri. The imaginary part of 2nd component gives us R^ = a. 
Thus we find out that k* = (2cos~^ a, 0). Then Eq (15^ become 



{2o? + + 1c?Rx = 



and one can verify that E = Ri indeed satisfying this equation. In summary, at k*, we have 
R2 = 0, R3 = a and E^ = Ri and the Hamiltonian becomes 



=+1 



V 






iRi 


a 


-iRi 





ia 


a 


—ia 






(38) 



The 3 eigenvalues are Ri and ± ^/Rj + Sa"^). Clearly, if i?i > and a is small 

positive number, Ri is top band. |(— -Ri + ^/ Rf + 8a^) ^ is the middle band. |(— -Ri — 
>y Rl + 8a^) ~ — -Ri is the lower band. So it is consistent with our previous assumption that 
Ri is top band. The eigenstate corresponding to E^ = Ri is 

\u+) = {t, 1, 0)^ (39) 

Now we can also compute the top eigenvalue of i?s^=_i at k*. At this point, we find 

/ n .0 . \ 



\ 



iRi a 
-iRi —ia 
a ia 



(40) 



The 3 eigenvalues from bottom to top are — -Ri, |(-Ri — -\/Rf~-i-"8a^)and |(-Ri + ^/R^~+~8a?) 
respectively. So the top band is E^ = |(-Ri + \/ Rf + 8a^) and the corresponding eigenstate 
is 

(41) 



\u.) 



\ 



2a^ + Rl + Riy/Rl + 8a^ 
-i(2a^ + Rj + Ri^Rl + Sa"^) 
a{3Ri + y^^TSa?) 



J 



Then we find ('U_|9|m+) oc j = 0. Thus ±k* is pair of zero points of -P(k). 

Now we show that this is the only pair of zeros of -P(k). Above we only discussed the 
special momentum point k* on which the wave function Eq (|37|) is not well defined. Other 
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than these points, we can always use Eq (137|) . Denote the top band of Hs^=i to be Ep and 
top band of Hs^=^i to be En- It is easy to find that 

P(k) = (M_(k)|e|M+(k)) oc (M_(k)|(|M+(k)))* = A + iB (42) 
A = a^{EpEn — R\) — Rl{Ep + E^) + E'^E'^ — EpEn{R\ + -R3) + R\R^ 
B = lai^EpEn — R^R^ 

Now we want to show A and B cannot be zero at the same time. 

First we consider the case that R2 7^ 0. If we assume 5 = 0, then this means EpE^ = R\- 
We know that Ep^n satisfy equations E^^^ — (a^ + R^)Ep^n ± 2aRiR^ = 0. Summing up 
these two equations and making use of the fact that Ep^n are positive, we find the following 
identity 

El + El = EpEn + a^ + R^ (43) 

Combine with EpEn = Ri, we find E^ + E^ = + R^ + R\. Now plug all the above 
result into the expression of A, we find that A = ~R\{o? + R\ + R^). We know that both 

+ Rl + Rl> and Rj = EpEn > 0, thus A ^ in this case. 

Second, we consider the case R2 = 0. From identity Eq. we find that 

A = {EpEnf + {a" -Rl- Rl){EpEn) - 2a^Rl (44) 

We also have Ep^n satisfying E^^^ — (a^ + R^) ± 2aRiR3/Ep^n = 0. Subtracting these two 
equations, we find that (Ep — En) = —2aRiR^ ^^^ Squaring the above equation and making 
use of Eq. ( H3|) , we find another identity 

o? + R^- EpEn = {2aRiR^f—^—. (45) 

Now if we assume that A = 0, then we find the following two equations 

{EpEnf + (a^ -Rl- Rl){EpEn) - 2a^Rl = (46) 
-{EpEnf + {a^ + R^){EpEnf - {2aRiR,f = (47) 

Multiplying Eq. ^ with {EpEn) and adding it to Eq. we find {EpEnf = 2RlRl + 

Rl{EpEn). Combining this resuh with Eq. (HHD, we find (a^ - Rl){EpEn - 2Rj) = 0. Since 
we do not consider the special point which satisfy R2 = and R3 = ±a, we find EpEn = . 
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By plugging this result back to Eq. (1461) . we will find a contradiction. Therefore in this case, 
we also have A ^ 0. 

In Summary, other than the special points, there is no more zeros of -P(k). Thus the 
Z2 invariance is —1 indicating the QSH state. Similar to the two band model, the singular 
points of the eigenfanuction and their gauge dependence can be explicitly demonstrated by 
connecting the Z2 invariant to the second order Chern number. But the calculation will be 
much more complicated than the two band case. 

VI. CONCLUSION 

We have discussed the the possible QSH state in a three-band model motivated by the 
lattice of copper-oxygen planes of the high-temperature superconductors. The calculation 
of the Chern parity or Z2 invariant is through a detailed analysis of the singular points of 
the eigenfunctions. We show that the the zeros of the PfafRans of time reversal operator 
are coincident with the singular points of the eigenfunction. It is known through the di- 
mensional reduction that the nontrivial Chern parity corresponds to the odd second order 
Chern number of a corresponding higher dimensional model. From a wave function point 
of view, we show that the geometric meaning of these topological invariant manifests in the 
obstruction which prevent us to define the wave function and Berry phase globally, just as 
the first order Chern number in the QAH state. An interesting fact is that the non-abelian 
Berry phase of this 4D hamiltonian can be mapped to a non-abelian gauge theory instanton. 
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